Power law dependence of the angular momentum transition fields in few-electron 

quantum dots 
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We show that the critical magnetic fields at which a few-electron quantum dot undergoes tran- 
sitions between successive values of its angular momentum (Af), for large M values follow a very 
simple power-law dependence on the effective inter-electron interaction strength. We obtain this 
power law analytically from a quasi-classical treatment and demonstrate its nearly-universal valid- 
ity by comparison with the results of exact diagonalization. 
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I. INTRODUCTION 

The field of quantum dots has been evolving under the 
busy attention of theorists and experimentalists alike for 
already almost 20 years by nowiSi^ This amount of 
interest and the enthusiasm is due to the unique blend 
of technological advances and fundamental interest. A 
special role is assigned to the magnetic field as a ver- 
satile tool with which one can tune the electronic prop- 
erties of the quantum dots. In particular, the applica- 
tion of a perpendicular magnetic field aids the Wigner 
crystallization;^ the formation of strongly correlated 
many-body states^S and induces ground state multiplic- 
ity transitionsiS 

From the theoretical side, the 'exact' numerical di- 
agonalization niethoc&i2iiiii24iiiiiiSii& is of profound 
importance as a reference for approximate treatments. 
With the presently available computing power this ap- 
proach has been successfully applied to compute the elec- 
tronic structure of few-electron quantum dots. However, 
this method also has sharp limitations since the numer- 
ical effort grows exponentially with the number of elec- 
trons. Therefore, the construction, refinement and appli- 
cation of less demanding approximate approaches based 
on some enlightening or visualizing idea is of critical im- 
portance. 

Among such ideas is the successful introduction of 
quasiparticles that came to be known as composite 
fermionsm^'i^ This transformation reduces the initial 
problem to a much more tractable problem of non- 
interacting or, at least, weakly interacting particles. An- 
other fruitful idea is to gain insight from the classical 
picture. The electrons in a quantum dot can be viewed 
as classical particles vibrating around their equilibrium 
positionsiSiSfliSi not unlike atoms in a molecule or a solid, 
while the molecule can rotate as a whole. These devel- 
opments culminated in the formulation of the rotating 
electron molecule approachi22i24 based on a classical visu- 
alization and successfully competing with the composite 
fermion model. 

In the present work, we call attention to yet another 



manifestation of the classical nature of quantum dots in 
strong magnetic fields. Namely, we concentrate on the 
critical magnetic fields, i. e. the field strengths at which 
the ground-state angular momentum of the dot switches 
between its two successive values. Here, we found that 
these critical field values show a strikingly simple de- 
pendence on the effective Coulomb coupling strength. 
We compare estimates extracted from an unsophisticated 
quasi-classical model with the exact-diagonalization re- 
sults thereby demonstrating the viability of such essen- 
tially classical paradigms and presenting one more argu- 
ment in their favor. 

The paper is organized as follows. In Sec. II we present 
the exact-diagonalization results indicating that the crit- 
ical magnetic fields obey a simple power law. In Sections 
HI and IV we introduce and solve a simple quasi-classical 
model, and compare its predictions to the exact results 
in Sec. V. The paper ends with a concluding Sec. VI. 



II. EXACT DIAGONALIZATION 

In Fig. n the well-known energy spectrum of two elec- 
trons in a parabolic do1iS (the Zeeman energy is not 
included) is plotted as a function of the perpendicular 
magnetic field strength B. We work with dimension- 
less variables so that the system energy is measured in 
units of hujQ with uq being the characteristic frequency 
of the confinement potential, while the magnetic field 
strength is expressed as the ratio 7 = ujc/i^o of the cy- 
clotron frequency tOc — eB/mc to the above-introduced 
confinement frequency. The relative importance of the 
electron-electron Coulomb interaction is characterized by 
the dimensionless coupling constant A = ao/as- Here 
flo = y^h/mujQ is the characteristic dimension of the 
parabolic quantum dot, and ub = eh/me^ is the Bohr 
radius. Fig. ^ corresponds to a relatively large (i. e. 
strongly interacting) dot of A = 5, and we include the 
magnetic fields up to 7 = 4 and angular momenta up to 
|M| — 5. We observe that for positive magnetic fields the 
preferred values of the angular momentum are negative, 
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however, for the sake of convenience we will henceforth 
use the symbol M for its absolute values and dispense 
with the "— " sign. 



8 





0/ 


y 




\ V 








- V \ 









































12 3 4 

Y 

FIG. 1: The spectrum of a strongly interacting (A = 5) two- 
electron parabolic dot as a function of perpendicular magnetic 
field. The lowest term of each angular momentum M = . . . 5 
is shown. The Zeeman energy is not included. At the critical 
fields 7Af the angular momentum switches from M to Af + 1. 

The most conspicuous feature of this spectrum is the 
presence of certain critical values 7m at which the ground 
state of the quantum dot changes the absolute value of its 
total angular momentum from M to M + 1. This feature 
is brought about by the electron-electron interaction and 
persists for any number of electrons in the dot. 

We performed exact numerical diagonalization for 
quantum dots containing 2, 3 and 4 electrons (see Ref. 
[l3i for the description of the procedure) and calculated 
the critical magnetic field values for various effective 
Coulomb interaction strengths. For the two-electron 
quantum dot we included A values up to 20, and for the 
case of three electrons in the dot the upper limit was set 
at A = 10. The case of four electrons is numerically more 
demanding and we restricted the calculation to values 
up to A = 2. We note that the effective electron-electron 
interaction strength can be easily tuned by varying the 
strength of the confinement potential, and consequently, 
the quantum dot size. In Figs.|5[0]these critical magnetic 
field values are shown in a double-logarithmic plot as a 
function of the dimensionless coupling constant A. We 
display the magnetic field range between In 7 = (that 
is, LOc — ^0) and ln7 = 3 (t^c ~ 20Li;o)- 

Notice that plotted in a a logarithmic scale these func- 
tions are nearly linear, besides the linearity is more pro- 
nounced for large momentum M and large magnetic field 
values. This linearity of the above curves 

In 7a/ = au + hu hi A (1) 

implies the simple power law dependence 

1M^\^''- (2) 




FIG. 2: Critical magnetic field strength values for two elec- 
trons. The full lines delimit regions characterized by different 
ground-state angular momenta starting with Af = at the 
lower left corner and increasing in steps of AAf = 1 at each 
consecutive boundary. Note the monotonic decrease of inter- 
line distances with increasing A and 7. 
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FIG. 3: The same as in Fig. H but now for three electrons. 
An additional phase boundary separating regions of angular 
momentum Af = and Af = 1 that shows a different behavior 
is not shown. Note that widths of some regions corresponding 
to more stable states are considerably larger. 

We expect that such a plain dependence is a consequence 
of the classical nature of the electron system in the quan- 
tum dot in the limit of a strong perpendicular magnetic 
field (see e. g. Ref. l25[l as we will show in the next section 
where a simple quasi-classical description is developed. 

We observe that the phase diagram corresponding to 
the two-electron quantum dot (Fig. [SJ is more regu- 
lar that the two others, namely, the nearly-straight-line 
boundaries appear to be parallel and the widths of the 
consecutive regions are decreasing monotonically. In con- 
trast, the plots pertaining to the cases of three (Fig. O 
and four (Fig. ^ electrons include some regions of con- 
siderably larger width which indicate an enhanced sta- 
bility of the corresponding ground states. In particular, 
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FIG. 4: The same as in Fig. |21 but now for four electrons. 
An additional phase boundary separating regions of angular 
momentum Af = and M = 2 that shows a different behavior 
is not included. The regions corresponding to more stable 
states are wider. 



such regions correspond to the maximum density droplet 
states of angular momentum M — 3 for three electrons 
and M — 6 for four electrons. 

Moreover, we see that at the lower right corner of Figs. 
13 and 0] (three and four electrons, respectively) corre- 
sponding to low magnetic fields and strong Coulomb cou- 
pling there is a slight deviation from coUinearity between 
neighboring lines. In the case of the four-electron dot the 
phase boundaries 72 and 73 even merge together indicat- 
ing that for InA > —0.3 there is a direct transition from 
the state with M — 2 into the state with M = 4. Simi- 
lar transitions where the angular momentum changes by 
AM = 2 were also found for higher values of angular 
momentum, namely, there are transitions 10 — > 12 and 
12 14. These properties of the four-electron quantum 
dot were already pointed out and discussed in Ref. 15. As 
we will see shortly, the quasi-classical model will overlook 
the presence of such irregularities. 

We note that in each of Figs.|31and0|an extra phase 
boundary is present at low magnetic fields which was 
omitted. These disregarded boundaries — see Fig. 2(a) of 
Ref. O and Fig. 3 of Ref. |13| for phase diagrams of three- 
and four-electron systems, respectively — display a dif- 
ferent trend. They are of a purely quantum-mechanical 
nature and result from a specific distribution of parti- 
cles among the low-energy levels of a quantum dot, i. e. 
they are effects that can not be reproduced in a classical 
model. 



III. QUASI-CLASSICAL DESCRIPTION 

The 2D iV-electron system confined by a parabolic 
quantum dot and placed in a perpendicular magnetic 
field B is described by the following dimensionless 



Schrodinger equation: 

{ff-£;}*==0, 

where 
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We consider the ground state energy of the state of 
fixed angular momentum M as a function of Coulomb 
coupling constant A and the magnetic field 7 



E = E{X,^,M). 
Solving the equation 

E{X,j,M) = E{X,j,M+l) 



(5) 



(6) 



for 7 will define the above-introduced critical magnetic 
field values corresponding to the angular momentum 
transitions from M to Af -|- 1. 

It is well known that in the strong magnetic field 
regimeS^ the electrons form a Wigner crystal which for a 
system up to five electrons is just a ring of equidistantly 
placed electrons. 

Due to the rotational symmetry around the z-axis the 
Hamiltonian Q commutes with the total angular mo- 
mentum operator 



M 



N 



d 



(7) 



Therefore, the simplest way to calculate the ground en- 
ergy of the system is to exclude the magnetic field from 
the Schrodinger equation Q by means of the transfor- 
mation 



N 



* ^ exp i ^ ipn/N ^, 



71=1 



with the ensuing scaling of coordinates 

-1/4 



r„(l-^7V4)" 



(8) 



(9) 



which enables to transform the initial problem into the 
equivalent problem 



{iJ- £:o(Ao,M)}* = 
with the Hamiltonian 

N N 

71<771 

without a magnetic field. 



Ao 



'-71 '- 7n I 



(10) 



(11) 
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The coupling constants and the eigenvalues of both 
problems are related as 



A = Ao{l + (7/2)2} 



1/4 



(12a) 

EiX,j,M) = £;o(Ao,M)Vl + (7/2)2-7^/2. 

(12b) 

Here, we again used our earlier convention whence the 
symbol M stands for the absolute value of the angular 
momentum. Note that the eigenvalue Eo{Xq,M) indeed 
depends only on the absolute value of the angular mo- 
mentum, and Eq. (|12b(l is written in accordance with our 
previous convention for M. Combining Eqs. p2|l and © 
we arrive at the following equation 



7j\/ 



2x/l + (7m/2)^ 



= Eq{Xo,M+1)-Eo{Xo,M) (13) 



which expresses the critical magnetic field values in terms 
of the reduced problem ((Tn|l . 



IV. DERIVATION OF THE POWER LAW 
DEPENDENCE 

Now we turn to the solution of the eigenvalue problem 
defined by Eq. H1U|) which involves two peculiar features. 
One is that we need to solve for a highly excited state 
with a large angular momentum M that corresponds to 
the ground state of the original problem Q , and this fact 
complicates the solution. Also, according to Eq. (|12a|l . in 
the strong magnetic field limit Eq. (|10() has to be solved 
for small values of the coupling constant Aq. Therefore, 
the solution must be attainable by means of some per- 
turbative technique. 

Neglecting the electron-electron interaction term in the 
Hamiltonian Hll|l we obtain the zero order equation. Due 
to the rotational symmetry the following radial equation 
for the zero-order one-electron fimction can be written: 



1 d d 



1 



Here, the symbol m stands for the single-electron angu- 
lar momentum. As this angular momentum is large the 
eigenvalue can be estimated by approximating the effec- 
tive potential 



Kff(r) 



m 
2^ 



1 



(15) 



by a parabolic potential in the vicinity of its minimum 
at 



|l/2 



(16) 



Here we disregard the first-derivative term in Eq. H14|) 
as its inclusion gives only negligible corrections in the 



TO ^ oo limit. This leads to the following estimate for 
the one-electron energy 



eo{m) w Voff(ro) = |m|. 



(17) 



and the total zero-order energy of all electrons in the dot 
becomes 



Eq = N\m\ = M. 



(18) 



Let us now estimate the characteristic energies of the 
two possible types of motion of the electrons: the radial 
vibrations and the motion of electrons along the ring. 

The energy of the radial vibrations can be obtained by 
expanding the effective potential at the equilibrium point 
To- The second derivative of the potential 



Kff(^o) 



(19) 



is of order unity and does not depend on the orbital mo- 
mentum m. This fact enables us to disregard these vi- 
brations as their energy is small in comparison to the 
zero order energy eoim). Moreover, these vibrations are 
the same for any orbital momentum M (they are not af- 
fected by the electron-electron interaction which is even 
smaller, i. e. X/vq <^ 1), and thus they can be neglected 
when solving Eq. H13|) . 

The energy of the longitudinal electron excitations 
along the ring is even smaller. Actually, in zero order ap- 
proximation the electron state under consideration is de- 
generate with respect to these excitations since the same 
total energy of the electronic ring H18(l can be constructed 
from various angular momentum distributions among in- 
dividual electrons. Thus, this angular motion is highly 
affected by the weak electron-electron interaction. The 
result of this effect is well-know: it leads to the Wigner 
crystallization of the electrons along the ring. 

The first order correction to the eigenvalue Eq (Aq, M) 
is obtained by including the electron-electron interaction 
term, but calculated classically for the electrons posi- 
tioned equidistantly on a ring of radius rg. This energy 
correction reads 



A£;(Ao,M) = ^fN 



(20) 



where /a? is the Coulomb energy of the equidistantly lo- 
cated electrons on the unit radius ring, in particular. 



h = 1/2, h = V3, .A = 1 + 2V2. 



(21) 



Then, combining Eqs. (|18|l and (|20|l we obtain the final 
expression for the energy to first order in powers of Aq 



E{Xo,M) 



/ N 

M + Xo^—Jn. 



(22) 



Inserting this expression into Eq. H13|l . solving it to the 
lowest accuracy in 7"^ powers, and taking the average of 
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-0.64 



-0.66 



-0.68 



-0.70 



-0.72 




FIG. 5: The power law coefficients bm as a function of the 
angular momentum M. The quasi-classical (QC) treatment 
predicts a universal value bu = —2/3 for any number of elec- 
trons shown by the solid line. Different symbols indicate the 
exact values of the coefficients for quantum dots containing 
2 to 4 electrons. The two-electron result displays a mono- 
tonic behavior and is spline-approximated by the dotted line 
to guide the eye. 



the difference of two inverse square roots we obtain the 
result 



(23) 



Finally, taking the logarithm of Eq. (|23|l and comparing 
to Eq. we extract the expressions for the coefficients 



bM 



ln(2M+l)-iln7V-^ln/w, 
2 

~3' 



(24) 



V. COMPARISON WITH THE EXACT 
NUMERICAL RESULTS 

The quasi-classical solution predicts a universal value 
for the power-law index 6m = —2/3 independent of the 
number of electrons. The comparison with the exact di- 
agonalization results is displayed in Fig. |31 As discussed 
above, in the four-electron dot we encountered two cases 
when a single AM = 2 transition occurs in place of ex- 
pected two AM = 1 transitions in the whole considered 
range of A and 7 values. In these cases (namely, the tran- 
sitions 10 12 and 12 — > 14) we used the M-value cor- 
responding to the average of the two expected AM — 1 
transitions. 

In general, the discrepancy with our analytical model 
is small as all results fall into a range between —0.64 and 
—0.72. As expected, the deviation of the exact result 
from the quasi-classical limit rapidly decreases with in- 
creasing angular momentum M , and already at Af = 8 
the difference is less than 2%. 



CD 




FIG. 6; The coefficient ajv/ as a function of the angular mo- 
mentum for the case of A'' = 2, 3 and 4 electrons in the dot. 
The full lines denote the quasi-classical (QC) result, and the 
symbols correspond to the results from exact-diagonalization. 



The exact results corresponding to the two-electron dot 
show a monotonous behavior as indicated by the dotted- 
line in Fig. [S] The reason is that the two-electron system 
in a parabolic dot possesses only one non-trivial degree 
of freedom as the centre-of-mass motion can be sepa- 
rated from the relative motion and is not excited in the 
ground state. In contrast, the internal motion of three- 
and four-electron dots involves more degrees of freedom 
and is more intricate. This results in a more complicated 
dependence of 6m on the angular momentum M . In par- 
ticular, enhanced non-monotonicitics are discernible in 
the vicinity of more stable maximum-density-droplet or 
"magic" states at M = 3 and 6 for three electrons and 
M = 6 for four electrons. 

Fig. El shows the comparison for the coefficient om- 
The quasi-classical model (|24ll predicts a logarithmic de- 
pendence of om on the angular momentum M . We see 
that for the two-electron system the correspondence is 
nearly perfect. 

The behavior concerning three and four electrons is 
more complex. At low values of the angular momentum 
(and thus low magnetic fields) there are notable discrep- 
ancies between the these sets of results in Fig. In 
particular, we observe wide "gaps" , i. e. abrupt jumps in 
the ttM values obtained from the exact-diagonalization 
at angular momenta M = 3 (M = 6) corresponding to 
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maximum density droplet states in three (four) electron 
systems. Nevertheless, the overall trend of these depen- 
dences still follows the quasi-classical prediction rather 
closely, and the discrepancy between the classical and 
fully quantum mechanical results vanishes very quickly 
with increasing M. 

VI. CONCLUSIONS 

In conclusion, we presented a quasi-classical theory of 
the magnetic properties of few-electron quantum dots 
whose main result is a simple power law dependence of 
the critical magnetic fields on the Coulomb coupling con- 
stant. While the main virtue of this theory lies in its 



relative simplicity and ability to provide results without 
resorting to heavy computation, the comparison of its 
predictions to the exact results also proves its robustness 
even in the realm of quantum mechanics. 
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